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We perform a global analysis of the experimental data of the electromagnetic nucleon form factors, 
in space-like and time-like regions. We give the expressions of the observables in annihilation 
processes, such as p -\- p ^ £~, t = e or yi,, m terms of form factors. We discuss some of 

the phenomenological models proposed in the literature for the space-like region, and consider their 
analytical continuation to the time-like region. After determining the parameters through a fit on 
the available data, we give predictions for the observables which will be experimentally accessible 
with large statistics, polarized annihilation reactions. 
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I. INTRODUCTION 

Elastic hadron electromagnetic form factors (EEs) are fundamental quantities for the understanding of nucleon 
structure. They contain information on the nucleon ground state, and constitute a further severe test for the models 
of nucleon structure, which already reproduce the static properties of the nucleon, such as masses and magnetic 
moments. Moreover, the dependence of FFs on the momentum transfer squared, should reflect the 

transition from the non perturbative regime, where effective degrees of freedom describe the nucleon structure, to the 
asymptotic region, where QCD applies. 

The magnetic proton FF, which is the dominant term in the elastic ep cross section, has been measured at values 
up to 31 GeV^ in the space-like (SL) region [1], and from pp or e+e“ annihilation up to 18 GeV^ in the time-like 
(TL) region [2]. Large progress has been recently done in the determination of the electric and magnetic proton form 
factors, based on the idea, firstly suggested in Ref. [3], to measure the polarization of the recoil proton in ep elastic 
scattering, when the electron is longitudinally polarized. 

Experiments, based on this method, have been performed at JLab up to = 5.6 GeV^ [5, 6]. A similar method, 
applied to the reaction d{e, e'n)p in quasi-elastic kinematics, has allowed the measurement of the neutron electric FF 
up to Q^=l GeV^ using a polarized deuteron target [7] and up to Q‘^=1A7 GeV^, measuring the polarization of the 
outgoing neutron [8]. The polarization method has been also successfully applied at low Q^, for a precise determination 
of the neutron FFs, at Mainz, and shows that Geu is definitely different from zero ([9] and refs therein). 

These results have been obtained thanks to the availability of high intensity, highly polarized electron beams and 
polarized targets, and to the optimization of hadron polarimeters in the GeV range. An extension of the measurement 
of the polarization transfer in e-|-p ^ e -|-p up to 9 GeV^ is in preparation [10]. 

More data are expected in future, in SL region, after the upgrade of Jlab, and in TL region, at Frascati and at the 
future FAIR facility at Darmstadt [4]. 

In the TL region [2], due to the poor statistics, the determination of FFs requires to integrate the differential cross 
section over a wide angular range. One typically assumes that the Ge contribution plays a minor role in the cross 
section at large and the experimental results are usually given in terms of |Gm|, under the hypothesis that Gb = 0 
or \Ge\ = |Gm|- The first hypothesis is an arbitrary one. The second hypothesis is strictly valid at threshold only, 
i.e., for r = q'^({Am?) = 1, but there is no theoretical argument which justifies its validity at any other momentum 
transfer, where q^ 4to^ (to is the nucleon mass). 

The measurement of the differential cross section for the process p + p ^ + i~ a.t a, fixed value of the total energy 

s, and for two different angles 9, allowing the separation of the two FFs, |GmP and |G_eP, is equivalent to the well 
known Rosenbluth separation for the elastic ep-scattering. However, in TL region, this procedure is simpler, as it 
requires to change only one kinematical variable, cos 9, whereas, in SL region it is necessary to change simultaneously 
two kinematical variables: the energy of the initial electron and the electron scattering angle, fixing the momentum 
transfer squared, Q^. Due to the limited statistics, the Rosenbluth separation of the \Ge\'^ and |GmP contributions 
has not yet been realized in TL region. Early attempts showed that the large error bars prevent to discriminate 
between the two hypothesis on |G_e| and |Gm| quoted above [11, 12]. 
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The \Gm\ values depend, in principle, on the kinematics where the measurement was performed and the angular 
range of integration. However, it turns out that these two assumptions for Ge lead to comparable values for [GmI- 

In the SL region the situation is different. The cross section for the elastic scattering of electrons on protons is 
sufficiently large to allow the measurements of the angular distribution and/or of polarization observables. Data on 
Gm are available up to the highest measured value, ~ 31 GeV^ [1] and this FF is often approximated according 
to a dipole behavior: 

Gm((3^)/Mp = Gd, with Gd=[l + , ml = 0.71 GeV^, (1) 

where Hp is the magnetic moment of the proton. 

It should be noted that the independent determination of both FFs, Gm and Ge, from the unpolarized e~ +p-cross 
section, has been done up to = 8.7 GeV^ [13], and the further extraction of Gm assumes Ge = Gm/^J- p■ The 
behavior of Gsp, deduced from polarization experiments, in which, more precisely, the ratio Gep/Gmp is directly 
related to the longitudinal and transversal component of the scattered proton polarization, differs from Gm/Mp, with 
a deviation up to 70% at Q^=5.6 GeV^ [6]. This is the maximum momentum at which new, precise data are available. 

The recent experimental data have inspired many new theoretical developments, and shown the necessity of a global 
representation of FFs in the full region of momentum transfer squared. 

FFs are analytical functions of q^, being real functions in the SL region (due to the hermiticity of the electromagnetic 
Hamiltonian) and complex functions in the TL region. The discussion of the constraints and consequences of a 
description in the full kinematical domain was firstly done in Ref. [14] and more recently in Refs. [15-18]. 

The extension of the nucleon models developed for the SL region to the TL region is straightforward for VMD 
inspired models, which may give a good description of all FFs in the whole kinematical region, after a htting procedure 
involving a certain number of parameters [19, 20]. 

The purpose of this paper is to update and compare some of the available models on the world data set in both TL 
and SL regions, and to predict time-like polarization observables, in framework of these models. The paper is organized 
as follows. In section H the expressions for the relevant polarization observables, in the process p + p £'^ + £~, 
£ = e or p, are given as a function of the electromagnetic FFs, in Section HI we update some of the fits of nucleon 
FFs on the available data, and discuss their extension to the TL region. In section IV we give the predictions of the 
considered models in TL region. 

II. OBSERVABLES IN TL REGION 

We develop a simple and transparent formalism for the study of polarization phenomena for p + p ^ £'^ + £~, in 
framework of one-photon mechanism. 

The calculations of the cross section and of the polarization observables for the process p + p^£~^+£~,£ = e 
or /i, in the annihilation channel are more conveniently performed in the center of mass system (CMS), Fig. 1. The 
momenta of the particles are indicated in the hgure and \ki\ = \k 2 \ = \k\ and \pi\ = \p 2 \ = \p\- Let us choose the z 
axis along the direction of the incoming antiproton, the y axis normal to the scattering plane, and the x axis to form a 

left-handed coordinate system. The components of the unity vectors are therefore p = (0,0,1) and k = (sin 6, 0, cos 9) 

with p-k = cos 9, where 9 is the electron production angle in CMS. The relevant kinematical variable is the antiproton 
energy, E, which is related to the four momentum transfer, = s = {ki + fe)^ = 4i?^, as, in CMS, ki + k 2 = 0. In 
the laboratory (Lab) system, one finds = 2m^ + 2mEL, where El is the Lab antiproton energy. The observables 
are calculated in the approximation of zero electron mass. 

The starting point of the analysis of the reaction p + p ^ + e— is the standard expression of the matrix element 

in framework of one-photon exchange mechanism: 

2 

M = ^ui-k 2 )'ypu{ki)u{p 2 ) F’iAr(g^)7M “-|^Tl 2 Ar(g^) u{-pi), (2) 

where pi, p 2 , ki and k 2 are the four-momenta of initial antiproton and proton and the final electron and positron 
respectively, > 4m^, q = ki + k 2 = Pi + P 2 - Pin and F 2 N are the Dirac and Pauli nucleon electromagnetic FFs, 
which are complex functions of the variable g^ - in the TL region of momentum transfer. 

In framework of one-photon exchange, the matrix element is written as the product of the leptonic and hadronic 
currents: 

2 2 2 

M = ^£pjf, = ^(4 Jo - £■ J) = J, 

q^ 


( 3 ) 
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FIG. 1: The kinematics of the process p + p ^ e + e”*" in the reaction CMS. 


p(p2) e*(k2) 



FIG. 2: One-photon mechanism for p -|-p ^ e + e'^ (with notation of four particle four-momenta). 


where ioJo = 0, due to the conservation of the leptonic and hadronic currents^. The expression for the leptonic 
current is: 


1 = 


1/^4 (o'- 


kd ■ k)(j)i, 


( 4 ) 


where pi(p 2 ) is the two-component spinor of the electron (positron), k is the unit vector along the final electron 
three-momentum, and for the hadronic current: 


J = \f^x\ 


GM{q^){d - pd ■ jP) + —j=GE{q^)pd ■ P 
\/t 


Xi, 


( 5 ) 


where xi and X’i are the two-component spinors of the antiproton and the proton, p is the unit vector along the three 
momentum of the antiproton in CMS. 


From this expression one can see the physical meaning of the particular relation between the nucleon electromagnetic 
FFs at threshold: 

GE{q^) = Gm( 4), 4 = 4m^. 

The structure pa • p describes the p-\- p annihilation from D-wave, i.e., with angular momentum 1=2. At threshold, 
where r ^ 1, the finite radius of the strong interaction allows only the S-state, and Gm(^^)- i=Ge{(1^) = 0. 


^ The conservation of the current implies that £ • g = 0, i.e., £090 — £ • g = 0, but g = fci + ^2 = 0 in CMS. Therefore, £090 = 0 for any 
energy go, be., £0 = 0. 
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From Eqs. (3), (4), and (5) one can find the formulas for the unpolarized cross section, the angular asymmetry and 
all the polarization observables. 


A. The cross section 


To calculate the cross section when all particles are unpolarized, one has to sum over the polarization of the final 
particles and to average over the polarization of initial particles: 


da 

dfl 


IMP k vm 


647r2g2 p’ 


m 2 


_ 1 g4 

\M\^ = -^^eabJab, Lb = Lit Jab = JaJ^ ■ 

Using the expressions (4) and (5), the formula for the cross section in CMS is: 


—i 

dnj^ 


(1 + cos^ 6 I)|Gm|^ + - sin^ dlGfil 


( 6 ) 


(7) 


where Af = — , a = e^/( 47 r) ~ 1/137, is a kinematical factor. This formula was firstly obtained in Ref. 

- 4m2) 

[ 21 ]. 

The angular dependence of the cross section, Eq. (7), results directly from the assumption of one-photon exchange, 
where the photon has spin 1 and the electromagnetic hadron interaction satisfies the P—invariance. Therefore, the 
measurement of the differential cross section at three angles (or more) would also allow to test the presence of 27 
exchange [ 22 ]. 

The electric and the magnetic FFs are weighted by different angular terms, in the cross section, Eq. (7). One can 
define an angular asymmetry, TZ, with respect to the differential cross section measured at 9 = ttI^, ctq [15]: 


/ da 

\d ^/0 


j = (To [1 -f TZcos^ 

/ 0 


where TZ can be expressed as a function of FFs: 


TZ = 


tIGm]^- [Gj;!^ 

tIGm]^ + \GeV 


( 8 ) 


(9) 


This observable should be very sensitive to the different underlying assumptions on FFs, therefore, a precise measure¬ 
ment of this quantity, which does not require polarized particles, would be very interesting. 

The dependence of the total cross section can be presented as follows: 


cr(g 2 ) = Af-TT 


21Gm1^ H— 


( 10 ) 


Polarization phenomena will be especially important 'm p + p ^ + (. . The dependence of the cross section on 

the polarizations Pi and P 2 of the colliding antiproton and proton can be written as follows: 


{Pl,P2) — f ^y{Ply p P2y) p AxxPixP2x p AyyPiyP2y + AzzPizP2z 


\dn 


PAxz{PixP2z P PizP2x)\^ 


( 11 ) 


where the coefficients Ai and {i,j = x, y, z), analyzing powers and correlation coefficients, depend on the nucleon 
FFs. Their explicit form is given in the following sections. The dependence (11) results from the P-invariance of 
hadron electrodynamics. 
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B. Single spin polarization observables 

In case of polarized antiproton beam with polarization Pi, the contribution to the cross section can be calculated 
as: 

Here the terms related to |G_eP and |GmP vanish. For the interference terms, the only non zero analyzing power is 
related to Pyi 

(S) ^hv = ^sm2eim{GMG*E). (13) 

When the target is polarized, one writes: 

Again the terms related to \Ge\'^ and |GmP vanish. Moreover, one can find A 2 = Ai = A. 

Eq. (13) has been proved also in Ref. [21]. One can see that this analyzing power, being T-odd, does not vanish 
in p + p ^ £+ +i- ', even in one-photon approximation, due to the fact FFs are complex in time-like region. This 
is a principal difference with elastic ep scattering. Let us note also that the assumption Ge = Gm implies Ay = 0, 
independently from any model taken for the calculation of FFs. 

C. Double spin polarization observables 

The contribution to the cross section, when both colliding particles are polarized is calculated through the following 
expression: 

^ Aab — mn 

where a and b = x,y,z refer to the a{b) component of the projectile (target) polarization. Among the nine possible 
terms, A^^y = Ay^ = A^y = Ay^ = 0, and the nonzero components are: 

{%) = sin2 0(^|GMp + i|GE|2jAf, 

= -sin"0(^|GMp-i|GsP^Af, 

(^) ^ (1-f cos^6>)|GmP - ^sin^ejGBp Af, 

{Si) A- = (s)/“ = ( 14 ) 

One can see that the double spin observables depend on the moduli squared of FFs, besides A,^^- Therefore, in order 
to determine the relative phase of FFs, in TL region, the interesting observables are Ay, and A^z which contain, 
respectively, the imaginary and the real part of the product GeG*^. 

III. RESULTS AND DISCUSSION 
A. The data 

The nucleon FFs world data were collected and listed in Table I for proton FFs and Table II for neutron FFs in SL 
region. 
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In Fig. 3 the nucleon FFs world data in SL region are shown: the ratio UpGep/Gmp (Fig. 3a) , the magnetic proton 
FF normalized to the dipole FF and divided by Hp (Fig. 3b), the electric neutron FF (Fig. 3c), and the magnetic 
neutron FF normalized to the dipole FF and divided by (Fig. 3d). For the electric proton FF, the discrepancy 
among the data measured with the Rosenbluth methods (stars) and the polarization method (solid squares) appears 
clearly in Fig. 3a. This problem has widely been discussed in the literature, (for a recent discussion see, for instance, 
Ref. [23]) and rises fundamental issues. If the trend indicated by polarization measurements is confirmed at higher 
[10], not only the electric and magnetic charge distribution in the nucleus are different and deviate, classically, 
from an exponential charge distribution, but also the electric FF has a zero and becomes eventually negative. This 
scenario will change our view on the nucleon structure and will favor VMD inspired models like [19, 20], which can 
reproduce such behavior. 

We included data issued from both kind of measurements in the fit, although if a consensus seems to appear that 
FFs extracted from polarization measurements are more reliable, as less affected by all kinds of radiative corrections. 
Our purpose here is not to get the best but to get a global description of the overall data. The precision and the 
number of points is very different for the different FFs, therefore one can obtain a good for a model that reproduces 
well, for example, the electric and magnetic FFs in the SL region and fails in giving the trend of Geu in TL region. 

We included in the fit the data on proton magnetic FFs which were published after 1973, and we did not include 
the data on the neutron electric FFs from Refs. [35, 36, 38, 41, 44] as data of much better precision were, later, 
available, in the same range. 

The data in the TL region are drawn in In Fig. 4a, b for the proton and in Fig. 4c, d for the neutron, respectively 
and summarized in Table III. As no separation has been done for electric and magnetic FFs, the data are extracted 
under the hypothesis that \Gen\ = \Gmn\- Concerning the neutron, the first and still unique measurement was done 
at Frascati, by the collaboration FENICE [45]. 


B. The models 

Among the existing models of nucleon FFs, we consider some parametrizations, which have an analytical expression 
that can be continued in TL region: predictions of pQCD, in a form generally used as simple fit to experimental 
data, a model based on vector meson dominance (VMD) [62], and a third model based on an extension of VMD, with 
additional terms in order to satisfy the asymptotic predictions of QCD [63], in the form called GKex(02L). We also 
considered the Hohler parametrization [64] and the Bosted empirical fit [65]. 

In order to help the reader, we report in the Appendix the explicit forms of the parametrizations previously 
published, with the parameters corresponding to the present fit, compared to the published ones. 

The pQCD prediction, based on counting rules, follows the dipole behavior (1) in SL region, and can be extended 
in TL region as [66]: 


[GmI 


m) 

In2(g2/A2)’ 


(15) 


where A = 0.3 GeV is the QCD scale parameter and A is a free parameter. This simple parametrization is taken to 
be the same for proton and neutron. The best fit ( Fig. 4, dashed line) is obtained with a parameter A{p)= 56.3 
GeV* for proton and A{n)= 77.15 GeV^ for neutron, which reflects the fact that in TL region, neutron FFs are larger 
than for proton. One should note that errors are also larger in TL region. 

A possible explanation of the fact that FFs are systematically larger in TL region than in SL region (which is true 
also in the proton case) is the presence of a resonance in the NN system, just below the NN threshold [67]. 

More pQCD inspired parametrizations exist for the form factor ratio Fb/Ei, which include logarithmic corrections, 
and have been recently discussed in Ref. [16]. However, some of these analytical forms have problems related to the 
asymptotic behavior. This will be discussed in a future paper. 

The analytical continuation to TL region of the other models is based on the following relations: 


Q = 




r ln((5^) = ln{q^) — in 

I 'JW = 


(16) 


Most of the models predict a different behavior for the electric and the magnetic FFs in TL region, whereas, as already 
mentioned, no individual determination of electric and magnetic FFs has been done yet. We chose to fit the data 
assuming that they correspond to the magnetic FFs for proton and neutron, Fig. 4a and 4c, respectively. Therefore, 
the curves for the electric FFs, in Figs. 4b and 4d have to be considered predictions from the models. Including or 
not the data on neutron FFs, in TL region, influence very little the fitting procedure. 
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The parametrization from Ref. [62] is shown as a dotted line, in Figs. 3 and 4. This model is based on a view of the 
nucleon as composed by an inner core with a small radius (described by a dipole term) surrounded by a meson cloud. 
While it reproduces very well the proton data in SL region (and particularly the polarization measurements), it fails 
in reproducing the large behaviour of the magnetic neutron FF in SL region. The present Ht constrained on the 
TL data and on the recent SL data does not improve the situation. In framework of this model a good global fit in SL 
region has been obtained with a modihcation including a phase in the common dipole term. However, the TL region is 
less well reproduced [20]. Therefore, the curves drawn in all the figures correspond to the original parameters, which 
give, in our opinion, a better representation of the whole set of data. 

The result from an update fit based on the parametrization GKex(02L) [63] is shown in Figs. 3 and 4 (solid line). 
It is possible to find a good overall parametrization, with parameters not far from those found in the original paper 
for the SL region only. The agreement is very good, for both proton and neutron FFs. 

The Hohler parametrization [64], contains also pole terms with adjustable parameters. The p-exchange contribution, 
however, is fully determined, with constants Hxed on -kN data. The model contains 17 parameters, already, so we 
did not try to readjust or refit the p-contribution. As noted in the original paper, such model is not suited to the 
extrapolation to TL region, because poles appear in the physical region. Constraining the parameters, in order to 
avoid these instabilities, worsens the description in the SL region. Therefore, we give only a Ht on all FFs, in SL 
region, Fig. 3 (dash-dotted line), corresponding to jndj ~ 1.7. The formulas as well as the original and updated 
parameters are also given in Appendix. Parametrization [19] can be considered a successful generalization, in TL 
region, based on unitarity and analyticity. It requires the modelization of ten resonances, five isoscalar and five 
isovector. 

The Hosted parametrization [65] is an empirical fit to nucleon FFs, in the SL region, based on simple formulas which 
are useful for fast estimations. It does not seem possible to find a unique function, which describes satisfactorily both 
the magnetic nucleon FFs and the electric proton FF, so the parameters are specific to each FFs. In the extension 
to TL region, as for the Hohler parametrization, one can not avoid poles and instabilities, and attempts to obtain a 
description in SL and TL regions remained unsuccessful. Therefore, we give the fit for the SL region, only, as dashed 
line in Fig. 3, and report in the Appendix the useful formulas and the updated parameters. As one can see from the 
table, they do not differ more than 20% from the published ones and the fit corresponds to /nd] ~ 2. 


IV. PREDICTIONS IN TL REGION 

We give the predictions for the cross section asymmetry and the polarization observables, for those models, described 
above, which give a good overall description of the available FFs data in SL and TL regions. The calculation is based 
on Eqs. (8), (13) and (14), for a fixed value of the angle d = 7r/4. 

As shown in Fig. 5, all these observables are, generally, quite large. The model [62] predicts the largest (absolute) 
value at ~ 15 GeV ^ for all observables, except A^z, which has two pronounced extrema. 

All observables manifest a different behavior, according to the different models. The sign, also, can be opposite 
for VMD inspired models and pQCD. The model [63] is somehow intermediate between the two representations, as it 
contains the asymptotic predictions of QGD (at the expenses of a large number of parameters). 

The fact that single spin observables in annihilation reactions are discriminative towards models, especially at 
threshold, was already pointed out in Ref. [69], for the process e’*' -I- e“ ^ p + p on the basis of two versions of a 
unitary VDM model. The present results, (Fig. 5), for the inverse reaction p + p ^ -\- e~ confirm this trend and 

show that experimental data will be extremely useful, particularly in the kinematical region around ~ 15 GeV 


V. SUMMARY 

The measurement of polarization observables and the possibility to access individual nucleon FFs in TL and SL 
regions at larger and/or with higher precision is foreseen in next future. 

A general analysis of the experimental data on nucleon electromagnetic FFs, extracted from elastic scattering and 
annihilation reactions, has been performed in the available kinematical region. 

Expressions of the experimental observables in the reaction p -\-p ^ -\- e~ have been derived in terms of the 

electromagnetic FFs, as a function of the momentum transfer squared. 

Some of the models on nucleon FFs have been reviewed, extended in TL region and used to give predictions on 
experimental observables which should be useful to plan future experiments. 

Many questions are still open. Recent data in the SL region show that the ratio Gep/Gmp deviates from the 
expected dipole behavior. In the TL region, the values of [Gm], obtained under the assumption that Ge = Gm, 
are larger than the corresponding SL values. This has been considered as a proof of the non applicability of the 



Phragmen-Lindelof theorem, (up to s=18 GeV^, at least) or as an evidence that the asymptotic regime is not reached 
[14]. The presence of a large relative phase of magnetic and electric proton FFs in the TL region, if experimentally 
proved at relatively large momentum transfer, can be considered a strong indication that these FFs have a different 
behavior. In particular, it will allow a test of the Phragmen-Lindelof theorem [14]. 

Large progress in view of a global interpretation of the nucleon FFs is expected from future experiments with 
antiproton beams: it will be possible, at the future FAIR facility at GSI, to separate the electric and magnetic FFs in 
a wide region of s and to extend the measurement of FFs up to the largest available energy, corresponding to s ~ 30 
GeV2. 

The angular distribution of the produced leptons will allow the separation of the electric and magnetic FFs. The 
measurement of the asymmetry TZ (from the angular dependence of the differential cross section for p + p ^ + £~) 

is sensitive to the relative value of |Gm| and jG^;]. In particular, the 0-dependence of the single spin and double spin 
polarization observables is very sensitive to existing models of the nucleon FFs, which reproduce equally well the data 
in SL region. 

Similar information can be obtained from the final polarization 'm -\- £~ p -\-p [69], but in this case one has to 
deal with the problem of hadron polarimetry, in conditions of very small cross sections. 

Only the study of the processes p + p ^ -\- £~ and p -\-p ^ + £~, [70, 71] will allow to measure 

proton FFs in the unphysical region (for s < Arri^, where the vector meson contribution plays an important role) and 
to determine the relative phase of pion and nucleon FFs. 


VI. APPENDIX 

The Sachs FFs are expressed in terms of the Pauli and Dirac FFs as: 

Gi = <(Q^) + rAf (g2), = <(g2) + 

One can introduce the isoscalar and isovector FFs F-® and F)’', i = 1, 2 as: 2F[ = F® -|- FI", 2F” = F® — F". 

Then, the isoscalar and isovector currents can be parametrized in terms of meson propagators, effective FFs, and/or 
terms which insure specific properties, according to the different models. 


A. Model from lachello, Jackson and Lande [62] and lachello and Wan [18] 


FFs are parametrized following the work [62] , with a modification that consists in adding a phase in the dipole 
term, g{Q^), for the extension in TL region. 


with g{Q^) = 


FUQ") 

F/(g2) 

F®(g") 
F^iQ^) 
1 


(1 -I- 7e*^g2)2 


gm 

2 

g(Q^) 

2 

gm 

2 

g{Q^) 


(1 ~ Plo — f3(p) + A 
(1 - Pp) + P, 






4 + 


Pp + SFp^^/tt 


(gp + Q^) + (4/i^ + g^)rpa(g2)//x^ 

2 


i \ 1 \ Foj I F(j) 

[Pp + Pn - I - _|_ Q2 + ,,2 _ 1 _ n2 


4 + Q" 


[Pp - Pn-^) 


pI + 8TpP-,,/tt 


{pl + Q^) + iApl + Q^)Tpam/p. 


and Q;(g^) = — 

TT 



vWT4/4) + y^ 


2^71- 


, with the standard values of 


the masses m = 0.939 GeV, pp = 0.77 GeV, p,^ = 0.78 GeV, p^j, = 1.02 GeV, = 0.139 GeV and the p width 
Fp = 0.112 GeV. 

The values of the six parameters are given in Table V. 
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B. Model from Lomon [63] 


^ [ 1.0317 + 0.0875(1 + QV0.3176)-^ , y mg, ] 

’ 2[ (l + g2/0.5496) Z+m^,+Q2j iWj + 


9p'\ uD(n,2\ 


1-1.1192- 


F^{Q^) = ^ 


N 5.7824 +0.3907(1+ QV0.1422)-i ^ 

2 (1 + g2/o.5362) ^ m2, + g2 2 W J 


^ .. 9p' \ rpD/r^2x 


6.1731--Kp.^jF2^(g2), 




9 ^ _”2+__ \ ^ _F<^/'n2 


f- 2^02 + ^-'7^ 2 r n2 ^2“(Q^) + 2 7n2 -^2^(Q ) + 

U mi + g2 m2, + g2 7 Umi + g2 


JUJ Juj' J4> J 


F7’'°(g^) = a = p, LV and Ai^d = Ai fori^", Ai,£, ee: A^, /or 

A2^^ + g2 Ai + g2 




The set of parameters is reported in Table IV. 


A? Q^ + pI 


g" = g' 


Zn[(A2, + g2)/A|^J 

ln{Aii/AQCD) 


C. Model from Hohler [64] 


This model is also based on a VMD parametrization: 


Ff(g2) = 0.5 0.955 + 


(l + g2/0.355)M l + g2/0.536’ 


F7(g2) = 0.5 5.335 + 


(l + g2/0.268)" l + g2/1.603’ 




+;(")(g2) = FZ(g2) + v^^^-. 

Y5f) + g2 


The parameters are given in Table VI. 
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D. Model from Bosted [65] 

The analytical expressions are inverse of polynomes as functions of Q, whereas Geu is described by a different 
function, as suggested by Galster [68]: 


= 




(17) 


G 


n _ 

E — 


a^JLnTGD{Q'^) 

1 + /3t 


(18) 


with al, a and /3 free parameters. In the present notation j = 1,2,3 corresponds to Gep and Gmu and Gmp, 
respectively. The inverse polynomes are of fourth order for Gep and Gmu and of fifth order for Gmp- 
The parameters are given in Table VII. 


a 




a 




FIG. 3: Nucleon Form Factors in Space-Like region: (a) proton electric FF, scaled by ppGMp (b) proton magnetic FF scaled 
by t-LpGo , (c) neutron electric FF, (d) neutron magnetic FF, scaled by finGo- The predictions of the models are drawn: from 
Ref. [62] (dotted line), from Ref. [63] (solid line), model from Ref. Ref. [64] (dash-dotted line), from [65] (dashed line). 
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FIG. 5: Angular asymmetry and polarization observables, according to Eqs. (13) and (14), for a fixed value of 0 = 45°. 
Notations as in Fig. 4. 
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Reaction 

Q" [GeV^] 

Observables 

Laboratory 

Year 

Reference 

e + p ^ e + p 

4.08-9.59 

Gmp 

DESY 

1966 

Albrecht et al. [24] 

V 

0.16-0.85 

G Ep : G Mp 

SLAG 

1966 

Janssens et al. [25] 

V 

0.69 - 25.03 

G Mp 

SLAC 

1968 

Coward et al. [26] 

V 

0.4 - 2 

Gsp) GMp 

Bonn 

1971 

Berger et al. [27] 

V 

0.670 -3.01 

G Ep : G Mp 

DESY 

1973 

Bartel et al. [28] 


0.999-25.03 

G Mp 

SLAC 

1973 

Kirk et al. [29] 

V 

2.86-31.2 

G Mp 

SLAC 

1993 

Sill et al. [30] 

V 

1- 3 

Gepi GMp 

SLAC 

1994 

Walker et al. [31] 

V 

1.75-8.83 

G Mp 

SLAC 

1994 

Andivahis et o7[13] 

V 

0.65-5.20 

Gepi GMp 

JLab, E94110 

2004 

Christy et al. [32] 

V 

0.65-5.20 

Gepi GMp 

JLab, Hall A 

2004 

Qattan et al. [34] 

e + p^e + p 

0.49-3.47 

Gep! Gmp 

JLAB, Hall A 

2000 

Jones et al. [5] 

e + p^e + p 

3.5-5.5 

Gep! Gmp 

JLAB, Hall A 

2002 

Gayou et al. [6] 


TABLE I: Data considered in the present analysis, for proton EEs, in SL region. 
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Reaction 

[GeV^] 

Observables 

Laboratory 

Year 

Reference 

ed —f epn 

0.04- 1.16 

Geti 

SLAG 

1965 

Hughes et al. [35] 

ed epn 

0.19,0.39,0.56 

Gmu 

New York 

1966 

Stein et al. [36] 

ed —> epn 

0.39-0.565 

Gmti, GEn 

DESY 

1969 

Bartel et al. [37] 

ed —> epn 

0.28-1.8 

Gmti 

Harvard 

1973 

Hanson et al. [38] 

ed epn 

2.5-10 

Gmti 

SLAG 

1982 

Rock et al. [39] 

^He 

0.16 

Geti 

MIT 

1991 

Jones-Woodward et al. [40] 

ed epn 

0.435-1.36 

Gmti 

New York 

1964 

Akerlof et al. [41] 

D{e, e'n)p 

.0.255 

Ge-ti 

MIT 

1994 

Eden et al. [42] 

D{e, e'n),D{e, e'p) 

0.125-0.605 

Geti 

Bonn 

1995 

Bruins et al. [44] 

D{e, e'n), D{e, e'p) 

0.235-0.784 

Ge-ti 

MAMI 

1998 

H. Anklin et al. [45] 

D{e, e'n)p 

0.15 

Ge-ti 

MAMI 

1999 

Herberg et al. ]47] 

D{e, e'n)p 

0.34 

G_b„ 

MAMI 

1999 

Ostrick et al. [48] 

D{e, e'n)p 

0.21 

Ge-ti 

NIKHEF 

1999 

Passchier et al. [49] 

^He{e, e'n)pp 

0.67 

Ge-ti 

MAMI 

2003 

Bermuth et al. [50] 

^Heie, e') 

0.1-0.4 

GEn, GMn , 

JLab 

2000 

Golak et aZ. [52] 

D{e, e'n)p 

0.495 

GEn 

JLab 

2001 

Zhu et al. [53] 

D{e, e'n)p 

0.5, 1 

GEn 

JLab 

2004 

Warren et al. [7] 

Die, e'n)p 

0.3-0.8 

GEn 

MAMI 

2004 

Glazier et al. [9] 

Die, e'n)p 

0.5-1.5 

GEn 

JLab 

1999 

Madey et al. [8] 


TABLE II: Data considered in the present analysis, for neutron FFs, in SL region. 


Reaction 

]GeV^] 

Laboratory 

Year 

Reference 

e"''e“ ^ pp 

4.3 

ADONE, Frascati 

1973 

Gastellano et al. ]54] 

pp e'^e~ 

3.52 

CERN 

1977 

Bassompierre et al. ]55] 

pp e''"e“ 

3.61 

CERN 

1983 

Bassompierre et al. ]56] 

e"'"e“ ^ pp 

3.75-4.56 

Orsay,DCI 

1979 

Delcourt et al. [57] 

e''^e~ —> pp 

4.0-5.0 

Orsay, DGI 

1983 

Bisello et al. [11] 

pp e'*"e“ 

8.9-13.0 

FERMILAB, E760 

1993 

Armstrong et al. [58] 

pp e'^e~ 

3.52-4.18 

GERN, LEAR 

1994 

Bardin et al. ]12] 

e"''e“ ^ pp 

3.69-5.95 

ADONE, FENIGE 

1994 

Antonelli et al. [59] 

1 

+ 

T 

8.84 - 18.40 

FERMILAB, E835 

1999 

Ambrogiani et al. [60] 

pp e''"e 

11.63- 18.22 

FERMILAB, E835 

2003 

Andreotti et al. [2] 

e'''e“ ^ nn 

3.61- 5.95 

ADONE, FENIGE 

1998 

Antonelli et al. [45] 


TABLE III: Data considered in the present analysis for TL region. 
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Parameters 

GKex(02L) 

This work 

9p' / fp' 

0.0608 

-0.0152 

Kp> 

5.3038 

-30.315 

Qui/ fui 

0.6896 

0.5994 

Kui 

-2.8585 

5.4557 

9<l>/f<i> 

-0.1852 

-0.287304 

«:</> 

13.0037 

18.0208 


0.6848 

0.802158 

g(.^')/fp' 

0.2346 

0.125397 

k{w') 

18.2284 

15.4868 

Ai 

0.9441 

0.776982 

A_d 

1.2350 

1.06593 

A 2 

2.8268 

3.64885 

^QCD 

0.150 

0.377101 

N 

1. 

0.84501 


TABLE IV: Parameters of the model from [63]. Fixed constants are = 3.706, its = —0.12, rrip = 0.776 GeV, rrisj = 
0.784 GeV, = 1.019 GeV, = 1.45 GeV, = 1.419 GeV. 
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Parameters 

[18, 62[ 

This work 

7 

0.250 

0.259 

Pp 

0.672 

0.757 

Puj 

1.102 

1.212 

P<k 

0.112 

-0.114 

OCcj) 

-0.052 

-0.028 

6 

53° 

47° 


TABLE V: Parameters for the model from Ref. [18, 62]. 


FFs 

par 

[64] 

This work 

FFs 

par 

[64] 

This work 

Fis 


0.61 

0.49 

Fiv 


1.49 

1.35 


b(‘>) 

0.94 

0.84 


bM 

^2 

4.33 

3.22 


^3 

3.20 

0.47 


^iv) 

^3 

8.47 

21.28 

Fis 

(l^s) 

ai 

0.75 

0.74 

Fiv 

(1,1^) 

0.06 

0.23 


(i.s) 

-0.61 

-0.57 


(1,4;) 

-0.32 

-0.55 


(i.s) 

“3 

-0.23 

-0.13 



0.10 

0.04 

F2s 


-0.15 

0.74 

F2v 


-2.06 

-1.93 



0.18 

-0.57 



0.23 

0.22 



-0.03 

-0.13 



0.23 

0.13 


TABLE VI: Parameters for the model from Ref. [64], corresponding to the fit in SL region, only. The starting set of parameters 
is from fit No. 9.1 


FF 

ap) 

[65] 

This work 

G Ep 

ai 

0.62 

-1.55 


Ct-2 

0.68 

10.40 


0-3 

2.80 

-11.53 


al 

0.83 

7.87 

Gmp! l^p 

of 

0.35 

0.28 


al 

2.44 

2.61 


ai 

1.04 

0.41 


a\ 

0.34 

1.03 


al 

0.34 

0.34 

Gmu j /^n 

Ol 

-1.74 

0.14 


al 

9.29 

3.39 


al 

-7.63 

-2.07 


al 

4.63 

3.09 

Geti 

a 

1.25 

0.43 


p 

18.3 

-0.22 


TABLE VII: Parameters for the model from Ref. [65], corresponding to the fit in SL region, only. 








